RKKY interaction in gapped or doped graphene 
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In our previous work (E. Kogan, Phys. Rev. B 84, 115119 (2011)) we calculated RKKY inter- 
action between two magnetic impurities in pristine graphene using the Green's functions (GF) in 
the coordinate - imaginary time representation. Now we show that the calculations of the GF in 
this representation can be simplified by using the Feynman's trick, which allows to easily calcu- 
late RKKY interaction in gapped graphene. We also present calculations of the RKKY interaction 
in gapped or doped graphene using the coordinate - imaginary frequency representation. Both 
representations, corresponding to calculation of the bubble diagram in Euclidean space, have an 
important advantage over those corresponding to calculation in Minkowskii space, which are very 
briefly reviewed in the Appendix to the present work. The former, in distinction to the latter, 
operate only with the convergent integrals from the start to the end of the calculation. 

PACS numbers: 75.30.Hx;75.10.Lp 



INTRODUCTION 

RKKY interaction between two magnetic impurities in 
graphene was theoretically studied quite intensely dur- 
ing last several years |ll4l9|. (A terse but precise review 
of the issue one can find in the book by M. Katsnelson 
[20l| . ) One may ask, why the problem, which is in princi- 
ple so simple (in the lowest order of perturbation theory, 
as it was treated in all the papers referenced above, the 
problem is equivalent to calculation of a single bubble 
diagram), was the subject of so many publication, using 
different approaches? One of the answers to this ques- 
tion is connected with the fact that a simply written inte- 
gral is not necessarily a simply calculated integral. More 
specifically, in most of the approaches mentioned above 
the integrals defining the RKKY interaction in graphene 
turned out to be divergent, and to obtain finite values 
from these divergent integrals one has to implement the 
complicated (and to some extent arbitrary) cut-off proce- 
dure. The cure from this disease turned out to be calcu- 
lation of this diagram not in Mikowskii, but in Euclidean 
space. Green's function (GF) in coordinate - imaginary 
frequency representation were used in Ref. 0. GF in 
coordinate - imaginary time representation were used in 
Refs. 0,[IH. 

In the present work we generalize the approach of our 
previous publication ll| to treat the case of gapped 
graphene. Not to distract attention of the reader from 
the aspects of the physics we are going to concentrate 
upon, we consider a toy model of graphene, with free 
electrons being described by the 2d Dirac Hamiltonian. 
The existence of two Dirac points in graphene leads to ad- 
ditional angular dependent factor in the formula for the 
RKKY interaction. This angular dependence was thor- 
oughly studied previously [j, 0] and does not interfere 
with the physics we are discussing in this work. 

The effective exchange RKKY interaction between the 
two magnetic impurities with the spins Si and S2, sitting 



on top of carbon atoms with the relative radius-vector 
R is 



(1) 



where J is the contact exchange interaction between each 
of spins and the graphene electrons, and x(-R) is the free 
electrons charge susceptibility, depending upon whether 
the carbon atoms belong to the same, or different sub- 
lattices. 



IMAGINARY TIME REPRESENTATION: 
GAPPED GRAPHENE 

The susceptibility expressed through the GF in the 
imaginary time -coordinate representation was written 
as flSEl 



x(R) 



-4/ 0(R;r)0(-R;-r)dr. 
Jo 



(2) 



Consider electrons in 2d described by Dirac equation 

[id t + ivj ■ d — 7°m] tp = 0, (3) 

where 7 • d = 7 1 <9i +J 2 d2- To find the GF in coordinate- 
time representation it is convenient to use the trick in- 
troduced by Feynman 2l|, that is to represent the GF 
as 



where 



G(R, t) = [id t +ivj-d + 7V] I(R, r), (4) 

T(R f) - [ d k 1 r -igfel*l+ik-R (K-) 



and Ek — V v 2 k 2 + to 2 . 

Going from real to imaginary time, instead of Eqs. 
and (O we obtain respectively 



G(R, t) = [d T + ivj-d + 7V] I(R, r), 



(6) 
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and 



I(R T ) = I — L e -S fc |r| +i k.R 



(7) 



Performing in Eq. ([7]) angular integration we obtain 



-in 



e -E k \r\ 



E, 



J (kR)dk 



(8) 



(Jo is the Bessel function of zero order). Using mathe- 
matical identity [22I ] 



00 ^—p^x' 2 +z 2 



we obtain 



\/x 2 + 



/(i?;r) 



Jo(cx)dx = 



z\/p 2 +£ 



( p2+c 2)l/2 



4tt« (uV'+.R 2 ) 1 ^ ■ 

In the gapless case mR/v C 1 we obtain [j| 
1 Tjr + 17 ■ R. 



a(ii;r) = - 



4tt (v 2 t 2 + i? 2 ) 3 / 2 ' 



(9) 



(10) 



(11) 



Taking into account that {7*, 7 J } = —2Sij (i = 1,2), and 
performing integration in Eq. @ we obtain [4| 



x cc (R) 



1 



x AB (i?) 



647rwi? 3 ' 



6AttvR 3 



(12) 



(CC can mean either AA or £?£>). 

In the opposite case mR/v > 1 we may approximate 
the GF as 



G(R;t) = 



m*y°R — «|m|7 • R 
AnvR 2 



-R\m\/v— vr 2 \m\/2R 



(13) 



Performing integration in Eq. ||2J), both for intra- 
sublattice, and for inter-sublattice susceptibility we get 



X 



1 / m \ 
~8v~ \ttvRJ 



.3/2 



-2R\m\/v 



(14) 



It is worth paying attention to the fact that Eq. (114)) 
seems to contradict rigorously proved theorem stating 
that for any bipartite lattice at half filling, the RKKY in- 
teraction is antiferromagnetic between impurities sitting 
on top of atoms belonging to opposite sublattices (i.e., A 
and B sublattices in graphene) , and is ferromagnetic be- 
tween impurities sittin g on top of atoms belonging to the 
same sublattice ji, 11, 23|. However, the theorem is not 
applicable to Hamiltonian ([3]), with its last term meaning 
that if we rewrite the Hamiltonian in the tight-binding 
representation, the intra-sublattice hopping will appear, 
hence the lattice is no longer bipartite. More specifically, 
the spectrum still has the symmetry of that in bipartite 
lattice, but the wave functions do not 11|. 



IMAGINARY FREQUENCY REPRESENTATION 

The approach will be based on equation Q 
1 f°° 

X (R) = — g 2 (R;u)dto, (15) 
n J -00 

where Q(R;uS) is the GF in the coordinate - imaginary 
frequency representation. 

Gapped graphene 



The Green's function is 

— ioj — v~/ ■ k — 7°m 
y{k,u3) = 5 5 5— T5 — . 

u z + m z + v z k z 

Hence for the diagonal part of the GF we obtain 



(16) 



Q {d) {RiUj) 



-lu-^m f°° kJ (kR)dk 



27T Jo w2 + m2 + V 2 k 2 



-ioj — 7 m 



-Kq (V^ 2 + m 2 R/v^j 



2irv 2 

and for the non-diagonal part 

-v f°° k 2 J 1 (kR)dk 



(17) 



g ab (r,oj) 



-Vw 2 + 



2irv 2 



27r J uj 2 + m 2 + v 2 k 2 

-Ki (V^ 2 +m 2 R/v S J . (18) 



Substituting into Eq. ([T5]) we obtain [2f 



x cc (R) = 



1 



47T 3 Z) 4 

(uj 2 - m 2 )Kl (yjtj 2 + m 2 R/v) du; 



(uj 2 + m 2 )K 2 (^uj 2 + m 2 R/v) du). (19) 

j 

For mR/v <C 1, using mathematical identity [13] 

poo 9 Q -3 /, 

/ x^K^K^dx = — — -r - 
' c a T(a) V 



a + ji + v 



ri^rnrf-rnrf -!!-'!, (20) 



we recover Eq. (| 1 2 p . 

For mi? /u» 1 we may use asymptotic expression for 
modified Bessel functions 



(21) 



After calculating the resulting integrals in Eq. (??) using 
the Laplace method, we recover Eq. (fl"4"|) . 
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Doped graphene 

For the case of doped (ungapped) graphene the GF is 

1 



(22) 



iui + fi — vj ■ k ' 

where fi is a chemical potential. From Eq. (|2"2"|) for the 
diagonal part of the GF we obtain 

— i(ui — ifi) f 00 kdk 



(27r) 2 J (uj - i/i) 2 + v 2 k 2 



271 



ikR cos 9 



d6 



— iji) 



kJ (kR)dk 



2n J (uj — i/i) 2 + v 2 k 

27rtr 

and for the non-diagonal part 

fc 2 dfc 



(23) 



<^ B (i?,c) = 



—v 



(27r) 2 y o (w - ill) 2 + v 2 k 2 



271 



ikR cos 6+i9 



d9= — 



k 2 J 1 (kR)dk 



27r J (uj — i[i) 2 + v 2 k 2 

sign(uj)(uj — ifi) 



lirv 2 



-Jfi[sign(w)(a;-i/i)E/u]. (24) 



(CC can mean either ^L4 or SB; i<Qi and i^i are the 
modified Bessel function of zero and first order respec- 
tively). We have used mathematical identity, valid for 
Re z > [12], 



o (a; 2 + z 2 )p 



J u (cx)dx 



c p-X z u-p+l 

2P~ 1 T(p) 



K u -p+i(cz). 



(25) 

The susceptibility p^]) is expressed through the inte- 
grals 



Re 



(w - vi)^o l I( w - in)R/v] duj \ . (26) 



Considering integrals in the complex plane it is conve- 
nient to deform the contour of integration and present 
the integrals as 

Re<! / uj 2 K 2 A (ujR/v)duj + / u> 2 Kl 1 (ujR/v)duj 



Taking into account the identity 

K a (-ix) = ?-i a+1 [J a (x)+iY a (a 



(27) 
(28) 



we get (12 1 

x° G {R) = x%{R) 



■k F R 

I - 1.6 / dzz 2 J (z)Y (z) 
o 



1G 



k F R 



dzz 2 Ji(z)Y 1 (z) 



(29) 



where kp = fi/v, and Xn=o{R) ar e given by Eq. (fT2"j) . 
The integrals in Eq. (1291) can be presented in terms of 
Meijer functions [l2|, [ijl (I address the reader to these 
References for the details). 

It is interesting to compare the RKKY exchange in 
doped graphene, with its two sublattices and linear dis- 
persion law, with that in ordinary two-dimensional elec- 
tron gas. For the latter the Green's function is 



Q(k,uj) 



1 



iui + ii — k 2 /2m 
and the susceptibility turns out to be [26[ 



1 



x(R) ~ W2 



k F R 



dzzJo(z)Yo(z). 



(30) 



(31) 



CONCLUSIONS 

In the end we would like to mention again that in the 
case considered, the GF calculations in Euclidean space, 
as it is not infrequently happens, have advantages over 
those in Minkowskii space. In particular, using the for- 
mer one have to operate only with the convergent inte- 
grals, in distinction to what happens when one uses the 
latter. 
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APPENDIX: APPROACHES BASED ON REAL 
FREQUENCY GF 

The approach, used in Refs. [l|, 0, 0] is based on equa- 
tion 



x(R) 



d 2 q 

(2tt) 2 



x(uj = 0,q)e 



where 



X,(uj = 0,q) 



d 2 k n F (g k ) - n F (6c+ q ) , 
(2tt) 2 £ k+q - E k ■■ 



(32) 



(33) 



= E n — H and uf(0 — (e^ + l) 1 is the Fermi dis- 
tribution function. This approach, though looking quite 
straightforward, brings with it a problem. In a model of 
infinite Dirac cones for \(ui = 0, q) we obtain a diverg- 
ing integral. To obtain finite values from these divergent 
integrals, as it was mentioned previously, one has to im- 
plement the complicated (and to some extent arbitrary) 
cut-off procedure Q ■ 
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The problem can be formulated in a different way. Be- 
ing calculated in a realistic band model, with the bands 
of finite width, x(cj = 0, q) is not a universal quantity. It 
depends not only on infrared physics, but on the proper- 
ties of electron spectrum and eigenfunctions in the whole 
Brillouin zone (even for small q). 

Another approach, formulated in Ref. [2{, starts from 
a well known equation for the susceptibility 



2i f°° 

X (R) = - G 2 (R, E)dE, 
n J-oo 



(34) 



where G is the retarded green's function. Here again the 
integral diverges on both limits of integration. However 
the authors changed the contour of integration, trans- 
forming the divergent integral (|34p into the convergent 
integral along the imaginary axis (see also Ref. [271]). 
The authors also considered RKKY interaction in gapped 
graphenc, when the power law decrease of the interaction 
with the distance turns into the exponential law. Actu- 
ally, the authors made the transition from real to imagi- 
nary frequencies, so some of our final results will be very 
close to ones obtained in Ref. [2} . 
The approach, using formula 



X{r,r') = 5n(r)/SV(r') 



(35) 



and, hence, calculating electron susceptibility on the ba- 
sis of equation 



x(R) 



Im [G 2 (R,E)] dE, 



(36) 



where Ep is the Fermi energy, was first used, in appli- 
cation to graphene to the best of our knowledge, in Ref. 
Q. An advantage of this approach is that it allows to 
easily consider the case of doped graphene, the disadvan- 
tage is that the approach, like the one presented above, 
has to deal with the divergent integral (the integral with 
respect to dE diverges at the lower limit of integration) . 
Also in this case, to obtain finite values from these di- 
vergent integrals one has to implement the complicated 
(and to some extent arbitrary) cut-off procedure. 
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